Close to the edge: hierarchy in a double braneworld 
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We show that the hierarchy between the Planck and the weak scales can follow from the tendency 
of gravitons and fermions to localize at different edges of a thick double wall embedded in an AdSs, 
spacetime without reflection symmetry. This double wall is a stable BPS thick-wall solution with two 
sub-walls located at its edges; fermions are coupled to the scalar field through Yukawa interactions, 
but the lack of reflection symmetry forces them to be localized in one of the sub- walls. We show that 
the graviton zero-mode wavefunction is suppressed in the fermion edge by an exponential function 
of the distance between the sub-walls, and that the massive modes decouple so that Newtonian 
gravity is recuperated. 
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I. INTRODUCTION 

The idea of confining our four-dimensional Universe 
inside a topological defect embedded in a higher dimen- 
sional spacetime dates at least as far back as the sugges- 
tion of Rubakov and Shaposhnikov [1] (see also [2] ) that 
we could be living inside a domain wall. They showed 
that fermions (of one chirality) can be confined to the wall 
by their Yukawa interactions with the scalar field. Do- 
main walls are vacuum configurations that are topolog- 
ically protected and thus stable, but their gravitational 
interactions cannot be ignored. This became evident af- 
ter the work of Randall and Sundrum [3] , in which they 
showed that the five-dimensional metric produced by an 
infinitely thin brane is enough to ensure that Newtonian 
gravity is reproduced on the brane. 

Actually, the RS brane does not even have to be a 
true domain wall, in the sense that no scalar field is 
needed for the effect. One just has to ensure that the 
bulk spacetime is Anti-De Sitter (AdS), and on the brane 
the four-dimensional cosmological constant is set to zero. 
The natural question is then whether fermions can also 
benefit from this effect, allowing us to live in a domain 
wall just because of its self-gravitation. The answer is 
on the negative, as shown in [4] fermion modes behave 
exactly opposite as the gravitons, the warp factor forc- 
ing them to escape from the wall into the bulk. If one 
wants matter to be confined to the wall, it is neces- 
sary to combine Rubakov-Shaposhnikov with Randall- 
Sundrum, i.e., to consider a real domain wall, made 
up of the vacuum expectation value of a scalar field 
that breaks a discrete symmetry, take into account its 
gravitational self-interactions, and make it couple to the 
fermions. This amounts to solve the five-dimensional 
coupled Einstein-Klein-Gordon system for an adequate 
potential, and many such solutions can be found in the 
literature. Among them, BPS walls, those where the 
four-dimensional cosmological constant is set to zero, 
are the most appealing "thick brane" generalizations of 
the Randall-Sundrum scenario. BPS solutions to the 
Einstein-Klein-Gordon system can be found by means of 



an auxiliary function of the scalar field, the fake super- 
potential of the first-order formalism of [5, 6, 7]. 

The fact that fermion zero mode localization requires 
exactly the inverse warp factor as graviton zero mode lo- 
calization, is used in this paper to provide a rationale for 
the large hierarchy between the Planck and weak scales 
in some particular thick wall solutions. These solutions 
are a straightforward generalization of the simplest BPS 
wall which represents a smoothing of the scenario of [3], 
found by Gremm [8]. The generalization produces an 
asymmetric double wall system: a BPS wall with a sub- 
structure consisting of two sub-wall located at its edges, 
but with different bulk cosmological constants on both 
sides. As a consequence of the asymmetry, fermions cou- 
pled to the scalar field are forced by the warp factor to 
localize on one of the sub-walls, the fermion sub-wall. 
On the other hand, the graviton zero-mode wavefunction 
is suppressed in the fermion sub-wall by an exponential 
function of the wall's thickness, i.e. the distance between 
the sub- walls, and gravitons localize in the opposite sub- 
wall, the Planck sub-wall. This allows one to provide a 
large hierarchy between the effective Planck masses on 
both sub-walls, as proposed by Randall and Sundrum in 
their earlier work [11] but with no orbifold geometries 
and no negative tension branes. 

Attempts to achieve suppressed mass scales with two 
positive tension branes, the so-called Lykken-Randall 
scenario [12], are found in the literature [13]. To our 
knowledge, however, all of them require some form of ra- 
dion field in order to stabilize the extra dimension [14]. 
In our case, the stability of the two wall system stems 
from their topological properties, they are just a special 
kind of BPS walls. The fact that no compactification of 
the bulk coordinate is required allows us to reproduce 
Newtonian gravity on the walls as in [3] , while keeping a 
large mass hierarchy as in [11]. Moreover, the fermions 
are not arbitrarily assumed to be located in a different 
wall as the gravitons, they do so as a consequence of 
spacetime being warped. In contrast with the scenario of 
[12], the fermion sub- wall is not a "probe" brane, and it 
has a non-negligible tension. 
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The paper is organized as follows. In the next Sec- 
tion we give an overview of the mechanism that provides 
stable, asymmetric double walls from a known BPS so- 
lution, and how fermions get localized on the brane with 
larger Planck mass. The following section is dedicated 
to explicitly construct these solutions from the simplest 
known BPS wall of [8], to find the graviton and fcrmion 
zero modes, and to calculate the Newtonian potential. 



II. DOUBLE WALLS AND HIERARCHY 

The so-called BPS double walls are solutions to the 
5-dimensional Einstein-scalar field set of equations that 
satisfy the BPS condition. As such, for the line element 
written in "proper length" coordinates 
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(1) 



they can be generated from a "fake superpotential" W{(j)) 
[5, 6, 7] by solving the BPS equations for the scalar field, 
warp factor and scalar potential 
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where prime denotes derivative with respect to the bulk 
coordinate ^. For W{(t>) given by 



a sm 



(3) 



with a and (t>Q real constants, solutions to (2) were found 
in [9] representing a family of double branes parametrized 
by an odd integer s > 1 that interpolate between AdS^ 
spacetimes. For s = 1, this is just the brane of Ref. [8], 
i.e. a regularization of the infinitely-thin RS brane, but 
for s > 1 the wall splits in two in a well defined sense: the 
energy density has two maxima as can be seen in Fig. 1. 
The fact that exponentiating the superpotential for a sin- 
gle wall gives rise to double systems was also used in [10] 
to construct BPS double walls. The topological charge is 
given by the asymptotic values of the superpotential 

= W{+oo) - W{-oo) = 2a (4) 

and is independent of s, i.e., it is the same as in the single 
wall. 

We have then a stable wall with two sub- walls, in which 
the thickness of the sub- walls goes to zero as s — > oo, 
while the separation between them remains fixed. The 
spacetime far away of the wall is AdS^ with the same cos- 
mological constant (A±) on both sides, while the space- 
time in between the sub-branes is nearly fiat, i.e. 
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A;. 
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The gravitational zero modes are calculated as usual 



with Ng a normalization factor. Because the spacetime 
between the two sub-branes is nearly fiat, the zero modes 
are not peaked at ^ = 0, but instead distribute smoothly 
over the whole system [15], as seen in Fig. 1. A similar 
behavior has been found for other BPS double branes 
[10]. 

Fermion modes of a given chirality can also be local- 
ized, by adding as usual a Yukawa coupling A with the 
scalar field [4]. One obtains 



4(C) 



-2A(5)„-A/0(C)d{ 



(7) 



The fermion zero modes for the system considered have 
been calculated in [16], and can be seen in Fig. 1. Details 
of the calculation are given in the next section, but the 
general behavior will suffice for the time being. 

Now, suppose the superpotential is shifted by a posi- 
tive constant 



(8) 



According to (2), we will have the same solution for the 
scalar field. Furthermore, since for the double wall (2,3), 
the extrema of V((t)) are the same ones of VF((/)), it follows 
that the extrema of the new scalar potential V'((/>) are 
the same ones of V^(</>), namely </> = and = ±7r0o/2. 
Hence, the two sub- walls are situated at the same place 
as before. Notice that the topological charge is the same. 
However, the cosmological constants arc now different at 
the two sides of the double wall and the spacetime in 
between the sub- walls is no longer flat, one has 

A_^-6(a-/3)^ k+^~&{a + l3f, Ajn -6/3^ 

(9) 

Accordingly, the new warp factor is from (2) 



A{0 = AiO - K 



(10) 



As long as we keep /? < a, the gravitational zero modes 
remain localized, but now they are exponentially sup- 
pressed for ^ > 



-2/3? 



(11) 



so that the graviton zero-mode function gets shifted to- 
wards the region with a smaller cosmological constant, 
i.e. with A_ . 

Fermion modes behave in exactly the opposite way, 
they are now 



4 e^K 



(12) 



so that they are shifted towards the larger cosmological 
constant region, i.e. with A+. In other words, the gravi- 
tons arc localized on the sub-branc situated at ^_ , hence- 
forth the Planck brane, while the fermion zero modes do 
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FIG. 1: Energy density (continuous), gravitational zero 
mode (dashed) and fermion zero mode (dotted) for the double 
branes with a > 0, /3 = 0, 5 = 4, s = 11, A = 4a, and arbitrary 
normalization. 



FIG. 2: Energy density (continuous), gravitational zero 
mode (dashed) and fermion zero mode (dotted) for the dou- 
ble branes with a > 0, 13 — 0.8a, S = 4,3 — 11, A = 4q, and 
arbitrary normalization. 



the opposite, they live in the opposite sub-brane around 
the weak brane. The situation is depicted in Fig. 2. 
The Planck masses in the Planck brane, M^l^'^'^'', and 
in the weak brane, Mp^i'^^^, are approximately related by 



M 



Planck 
PI 



(13) 



with A the inter-sub-brane separation. Since A here is a 
coordinate-dependent quantity, Eq.(13) should be made 
precise by calculating and comparing the gravitational 
potential on both sub-branes, which we do in the next 
section. 

As the thickness of the sub-branes approaches zero, 
this solution amounts to having a thick domain wall with 
two infinitely-thin sub-branes, the Planck brane and the 
weak brane, situated at different edges (in the transverse 
direction) of the wall. Clearly, the crucial role here is 
played by the double-brane solution, for which BPS sta- 
bility arguments apply [6]. 



III. EXPLICIT SOLUTIONS 



A. Double branes 



Solutions to the system (2) with un-shifted superpo- 
tential (3) (i.e. reflection symmetry) have been found 



in [9] by switching to the so-called "gauge" coordinates, 
where the metric is written as 



(14) 



The field equations can be integrated to give 
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cos2(0/0o) -2 
(15) 
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Hiy) 

V{cf>) = 3a^[sin((/)/0o)] 

, , fay 
<P[y) = (Poarctan^^ — 

The parameter S gives the thickness of the wall. For 
3 = 1, this is just the brane of Ref. [8] and it can be 
rigorously shown that the distributional limit (5^0 gives 
the RS spacetime [17]. But for non-zero 6 and s > 1, the 
energy density has two maxima, as can be seen in Fig. 1. 
Other double wall systems have been found in [10] and 
[18]. 

These walls have been studied in detail in a series of 
papers, their thin wall limit (5 — > in [9], localization of 
gravity in [15] and of chiral fermion modes in [16]. For 
6—1, there is a normalizable zero mode and a continu- 
ous of massive states that asymptote to plane waves as 
y — > ±oo. This continuum of modes decouple, in the 
sense that they generate small corrections to the New- 
tonian potential of the wall, at scales larger than the 
fundamental length-scale of the system, a~^. 

Let us consider the double domain wall (14, 15) for 
s ^ 1. The separation between the sub-branes, i.e. the 
distance between the maxima of the energy density, is 
given by 



1 



26 



l/2s 



(16) 



which behaves as 2^/a for s ;2> 1, while their thickness 7 
is approximately given by 



7---ln(|)+0(.-^) 
as Vo / 



(17) 



from which follows that 7 ^ for s :§> 1. Note that 
the double wall is lost in the (5^0 limit [9] , as expected 
from the fact that the distributional geometry 6 —>■ may 
be identified with the asymptotic behavior (i.e. far away 
from the wall) of the domain wall spacetime [19]. 

Following [17], the distributional limit s — > 00 of the 
Einstein tensor is found to be 

lim iGl + Ag%) = 

-3a [6iy - y+) + 6iy - y.)] {d'^dh + d^dxl) (18) 
with A given by 



A 
6" 



a^, \y\ > 6/a 
0, \y\ < 6/a 



(19) 
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Hence for s ^ 1 we have two infinitely thin walls located 
a.t y+ = 5 / a and = —5 /a, with tensions (t+ = ct_ = 
3a > 0. On the other hand, we find that the limit s ^ oo 
of the metric tensor gives the line element 



)26 / \ 2 

dy^ Ay\>6/c 



-dt^ + dx^dx' + dy^ 



, \y\ < 6/a 
(20) 

and the distributional Einstein tensor (18) turns out to 
be the Einstein tensor of the limit metric (20). 



B. Asymmetric double branes 

Now, shifting the superpotential as in (8) results in a 
solution with 



Hiy) 
My) 

T/(0) 

4>{y) 



H{y) 

A{y) - PyTiy) 

F(0) - 6/3[/3 + 2a(sin(0/0o))'-'/l 



(21) 



where Jf'(y) = 2-Fi[l/2.s, l/2s, 1 + l/2s, -(a?;/^)^'*] is the 
hypcrgcometric function. For |Q!y/(5| < 1, i.e. inside 
the double brane, and large s, 2F1 ~ 1 is a very good 
approximation, although we have used the exact solution 
whenever performing numerical integrations. 

The equation for the gravitational zero modes is found 
as usual. In the axial gauge h^y = 0, writing the 
transverse-traceless part of the metric perturbations /i^j^ 
as 



(22) 



one obtains 



l^^U{y)]i;,,^m'e-^^+^"i>,, (23) 



with 



U{y) 



-2A+H /2 f^^2A-H /2yi 



(24) 



As usual, in order to get a Schrodinger-like equation, one 
must change to conformal coordinates. However, from 
(23) we see that the zero mode is normalizable 



^2A-H/2 



(25) 



and since C/ ^ at |y| ^ 00, there is no gap between 
the massless and the massive modes. This can be seen 
by writing U in terms of W 



U{y) 



H/2i^^-H/2yi _ ^2H 



dW 

~d4 



(26) 



Clearly, shifting W hy & constant does not change the 
asymptotic behavior of U and the massive modes are 



therefore expected to decouple [20]. The zero mode, in 
proper length coordinates, is shown in Fig. 2. It follows 
that the graviton zero mode gets localized on the sub- 
wall located at the edge of the wall close to the region 
of smaller AdS^ curvature, thus defining the Planck sub- 
wall. A similar localization of the graviton zero mode 
on a rather different asymmetric double wall system was 
found in [18]. 



C. Adding fermions 

Fermion confinement in the double walls (14,15) was 
studied in [16]. In these coordinates, 5-dimensional 
spinors coupled to the scalar field by a Yukawa term of 
the form 



give confined chiral fermion modes on the wall 



-2l(y) -A J 0(y)e-f 



Uy 



(27) 



(28) 



for sufficiently large A. It was found in [16] that in general 
thin walls require large Yukawa couplings, which in the 
asymmetric case are bounded from below by the largest 
cosmological constant, in order to confine fermions. For 
the double-brane system (21), however, the Yukawa cou- 
pling depends also on 5 and we get 



A > 



2^ 



37r 



- 1/2) 



(29) 



Therefore, the Yukawa coupling can be kept at reasonable 
values if the thickness of the sub-branes 7 is decreased 
(recall that 7 — > for s ^ 1) while the separation be- 
tween them (the brane thickness ^ 5/a) is increased. 

The equation for the fermion zero modes was inte- 
grated numerically, and results are given in Fig. 2. As 
argued above, the fermion zero modes get located in the 
opposite sub-wall to the Planck sub-wall, thus defining 
the weak brane. 



D. Newtonian potential 

In order to have an effective four dimensional gravity 
on the weak sub-wall, we should demonstrate that the 
massive graviton modes decouple. Before presenting the 
results for the double wall system, it is instructive to 
consider the case s = 1, a single wall without reflection 
symmetry, with different cosmological constants A+ and 
A_ on both sides. A similar case was studied in [15] in 
the limit of a very small asymmetry. In the general case, 
the massive modes in the small m/ yJ7\Z approximation 
arc given by 



7,5/2 



A 



5/4 




(30) 
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with K of order one for any value of the parameters. In 
the reflection-symmetric case A_ = A+, this modes give 
the weU-known contributions to the Newtonian potential 
between masses separated by a 4-D distance r propor- 
tional to However, in the asymmetric brane the 
modes behave as m^/^, and their contributions are pro- 
portional to l/r*", the equivalent of having 6 extra com- 
pact dimensions. We expect that the double asymmetric 
walls exhibit a similar behavior. 

To calculate the contribution of the massless and mas- 
sive modes to the Newtonian potential, we switch to con- 
formal coordinates 



-dr + dx^dx' + dz^ 



(31) 



and consider the s ^ oo limit, i.e. the infinitely thin 
sub- wall idealization of (21). This limit is not as straight- 
forward as (20) for the symmetric case (15), due to the 
presence of the hypergeometric function in the warp fac- 
tor. However, it can be approximated by 



r [e+W" + (a + /3)(z-z+)] \ z+ 
3^'-''^ ^ {[cosh pS/a + 13 z]~\ 



< Z 

Z- < z < z+ 
{a — (3){z — z-)\ , z < z_ 

(32) 

with z± — ±13 ^sinh/35/a, < /? < a and 5 > 0. For 
/3 ^ we obtain (20), the limit s oo of the symmetric 
double wall (14, 15), written in conformal coordinates. 
The Einstein tensor of (31, 32) is given by 

G\ + Kg\ = -{d'^dtb + dl^dxl) 

[a+5{z-z+)+(T-5{z-z-)] (33) 

with A given in each region by (9) and where the sub- 
brane tensions are <7+ = 'iae^^^^" and cr_ = 3ae^'^''/". 
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Notice that although this resembles the two positive ten- 
sion three-branes scenario of [12], in our case the two 
branes separate AdS^, slices with different cosmological 
constants. Furthermore, the weak brane is not a so- 
called "probe" brane since its tension is not small (in 
fact, in proper length coordinates both sub- wall tensions 
are equal). 

Now, writing the transverse-traceless part of the metric 
perturbations h^^ as 



V(^,^) = e^^-"e^^^^/'<.(^) 



(34) 



with h^j^z ~ 0, one finds that the gravitational modes ijj^^, 
satisfy the Schrodinger equation 

(-^ + ^QM)^^.(^) = mV^,(z) (35) 



where 



Vqm ^ [e'^^'^/'Y e-^^<-'^/^ 



(36) 



For m = the solution is 
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sinh(2/M/a) 



The massive modes of (35) are given by 



-1/2 



(38) 



' (fc;i + \z~ z+|)i/2 [y2(m(fc;i + |z - z+D) + C+ J2(m(fc;^ + \z- z+\))] , z+ < z 

<(z) = N,i I (fcQ-i + z)i/2 [AY2imiko' +z)) + B J2{'m{k^^ + z))] , z. < z < z+ 

XkZ^ + \z~ z_|)i/2 [Y2{m{kZ^ + \z~ z_|)) + C-J2{m{kZ^ + |z - z_|))] , z < z_ 

I 



(39) 



where Y2 and J2 are the Bessel functions of order two, 
k± ~ (a ± /3) exp{=F/3i5/a}, fco = (3[cosh{(36 / a)\~^ and 
C-I-, C_, A, B are constants determined by the matching 
conditions. 

The contribution of the massive modes to the New- 
tonian potential can be now calculated by expanding 
around masses smaller than the smaller scale of the sys- 
tem, which is of order (3e~^^/°' (notice that the symmetric 
case /? = has to be treated separately). The normal- 
ized wave functions for the massive modes in the weak 



and Planck branes are found to be 



5 I a2 



1/2 



1/2 



-pS/a 



5/2 



(40) 



5/2 



^3/3<5/2q 



(41) 

to leading order in e , were k± are functions of a and 
(3 of order unity. Here, in order to obtain the correct nor- 
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malization constant for the wave function of the massive 
modes, we have used two regulator branes with positions 
at ±Zr with Zr > \z±\ and taking the limit Zr oo at 
the end of the calculation, extending to our two brane 
system, as close as possible, the single brane treatment 
of [21]. Wc can see that the modes behave as m^/^, as 
in the single asymmetric brane, and therefore for r < Tc 
the Newtonian potential behaves as in a scenario with six 
extra dimensions. 

The Newtonian potential in the weak brane is 



TOim2 



V+ir) 

r 

and in the Planck brane 

G^mim2 
V-{r) = 



1 



1 



r / 



(42) 



(43) 



to leading order in e~^^^", where Vc = exp {5(3S/3a) 
and the four-dimensional gravitational constants are 
given by 

^G^iN^fe-'P'/^, Gl^G^iNlfe^P"" (44) 
The ratio of the Planck masses on both branes is then 

M^'^^V^^p^""'' = e^'^^'" (45) 

and we can now refine our result (13). For Mpj''"'^'^ of 
order TeV, we would get a Planck mass in the weak brane 
of order 10^^ GeV by setting 



^^.12 



a 



IV. SUMMARY AND OUTLOOK 



(46) 



We have shown that the large hierarchy between the 
Planck and the weak scales can be attributed to the ten- 
dency of gravitons and fermions to localize at different 
edges of an asymmetric double domain wall. The embed- 
ding in a five dimensional Anti-de Sitter geometry which 
is not reflection symmetric determines which sub-brane 
is the weak brane and which one is the Planck brane. 



Fermions of one chirality are localized by their Yukawa 
interactions with the scalar field, but the warped metric 
forces them to live in a different wall than the gravitons, 
without need for additional assumptions. By calculat- 
ing the massive mode contributions in this system, we 
have shown that Newtonian gravity is recuperated in the 
sub-wall where the matter is located. 

The double-wall systems are straightforward general- 
izations of single, reflection symmetric BPS walls, and 
while here we have considered the case of the simplest 
one, the confinement of graviton and fermion zero modes 
on different sub- walls of the system is presumably generic 
to other double solutions without reflection symmetry. 
The generalization consists simply in taking a power of 
the fake superpotential, and then adding a constant to 
it. Since the double asymmetric walls are also BPS and 
therefore have a topological charge (which is in fact the 
same as the original wall), stability is guaranteed, and no 
additional fields or stabilization mechanism are required. 

While we believe this to be an interesting effect, many 
questions would have to be answered before attempting 
to use it as a solution to the hierarchy problem. For ex- 
ample, a mechanism for confinement of the gauge fields 
would be required, in particular, one that makes use of 
the fcrmionic fields to achieve localization such as in [22] 
would be well-suited, since then gauge fields will be con- 
fined to the matter wall. Scalar fields are also a problem, 
since the asymmetry drives them to the graviton's wall, 
their modes being proportional to the warp factor. An 
adequate coupling with the wall's scalar field could help. 
We hope to address this issues in a future publication. 
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